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ABSTRACT 



We compare the central mass concentration of Cold Dark Matter halos found 
in cosmological N-body simulations with constraints derived from the Milky Way 
disk dynamics and from the Tully-Fisher relation. For currently favored values of 
the cosmological parameters (Oq ~ 0.3; Aq = 1 — ~ 0.7; h ~ 0.7; COBE- and 
cluster abundance-normalized erg; Big-Bang nucleosynthesis O^), we find that halos 
with circular velocities comparable to the rotation speed of the Galaxy have typically 
three times more dark matter inside the solar circle than inferred from observations of 
Galactic dynamics. Such high central concentrations of dark matter on the scale of 
galaxy disks also imply that stellar mass-to-light ratios much lower than expected from 
population synthesis models must be assumed in order to reproduce the zero-point of 
the TuUy-Fishcr relation. Indeed, even under the extreme assumption that all baryons 
in a dark halo are turned into stars, disks with conventional /-band stellar mass-to-light 
ratios {M/Lj ~ 2 it 1{M/Li)q) are about two magnitudes fainter than observed at a 
given rotation speed. We examine several modifications to the ACDM model that may 
account for these discrepancies and conclude that agreement can only be accomplished 
at the expense of renouncing other major successes of the model. Reproducing the 
observed properties of disk galaxies thus appears to demand substantial revision to the 
currently most successful model of structure formation. 



Subject headings: cosmology: theory - galaxies: formation, evolution - methods: 
numerical 
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1. Introduction 

Over the past few years, cosmological models based on the paradigm of an inflationary 
universe dominated by cold dark matter (CDM) have proved remarkably successful at explaining 
the origin and evolution of structure in the Universe. The free parameters of astrophysical 
relevance in this modeling are surprisingly few: the current rate of universal expansion, Hq; the 
mass density parameter, Qq; the primordial baryon abundance, fib', and the overall normalization 
of the power spectrum of initial density fluctuations, as- Over the past few years, limits on the 
values allowed for these parameters have been consistently refined by improved observational 
techniques and theoretical insight, and it is widely accepted that a new "standard" model has 
emerged as the clear front-runner amongst competing models of structure formation. 

This model, which we shall call "standard" ACDM, or sACDM for short, envisions an 
eternally expanding universe with the following properties (Bahcall et al 1999): (i) matter makes 
up at present less than about a third of the critical density for closure (Qq 0.3); (ii) a non-zero 
cosmological constant restores the flat geometry predicted by most inflationary models of the 
early universe (Aq = 1 — ~ 0-7); (iii) the present rate of universal expansion is Hq ~ 70 km s~'^ 
Mpc~^ {h = Hq/IOO km s"^ Mpc~^ ~ 0.7); (iv) baryons make up a very small fraction of the mass 
of the universe (Oj, w 0.0125 ~ 0.0255 <C ^o)'-, and (v) the present-day rms mass fluctuations 
on spheres of radius 8 Mpc is of order unity (erg 1.1). The sACDM model is consistent 
with an impressive array of well-established fundamental observations: the age of the universe 
as measured from the oldest stars (e.g., Chaboyer et al 1998), the extragalactic distance scale as 
measured by distant Cepheids (e.g., Madore et al 1998); the primordial abundance of the light 
elements (e.g., Schramm & Turner 1998), the baryonic mass fraction of galaxy clusters (e.g.. White 
et al 1993), the amplitude of the Cosmic Microwave Background fluctuations measured by COBE 
(e.g., Lawrence, Scott &; White 1999), the present-day abundance of massive galaxy clusters (e.g.. 
Eke, Cole & Prenk 1996), the shape and amplitude of galaxy clustering patterns (e.g., Wu, Lahav 
Sz Rees 1999), the magnitude of large-scale coherent motions of galaxy systems (e.g., Strauss & 
Willick 1995, Zaroubi et al 1997), and the world geometry inferred from observations of distant 
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type la supernovae (e.g., Perlmutter et al 1999, Garnavich et al 1998), among others. 

Because its major parameters are fixed, sACDM is an eminently falsifiable model whose 
predictive power may be used to ascertain its validity on scales different from those used to tune 
the model. One scale of particular interest is that of individual galaxies, since few observational 
constraints on these small scales have been used to adjust the parameters of sACDM. This exercise 
is especially compelling because the dark halo structure found in cosmological N-body simulations 
of Cold Dark Matter universes seems at odds with dynamical studies of disk galaxies that assign 
a substantial gravitational role to the disk component (see, e.g., the "maximal disk" solutions of 
Debattista & Sellwood 1998, and references therein) as well as with rotation curve studies of dark 
matter-dominated galaxies (Moore 1994, Flores Sz Primack 1994, McGaugh Sz De Block 1998, 
Moore ct al 1999, Navarro 1999). These claims arc based largely on comparisons of the detailed 
shape of the rotation curve of very low surface brightness dwarfs with the innermost density 
profile of simulated dark halos. Unfortunately, the scales where deviations are most pronounced 
(the inner few kpc) are also the most compromised by numerical uncertainties (most simulations 
relevant to this problem published to date have gravitational softening scales of order 1-2 kpc). 
The comparison is thus rather uncertain. For example, Kravtsov et al (1998) have argued, on 
the basis of simulations similar to those used by the other authors, that CDM halos are actually 
consistent with the rotation curves of dark matter-dominated disks, a somewhat surprising result 
that illustrates well, nonetheless, the vulnerability of numerical techniques on scales close to 
the numerical resolution of the simulations. Prom a strictly pragmatic numerical standpoint, it 
would be desirable to circumvent these uncertainties by adopting comparison criteria that are less 
sensitive to the effects of numerical shortcomings. 

One possible choice is to use, rather than the dark matter density profile near the center, 
the total amount of dark mass within the main body of individual galaxies. For spiral galaxies, 
this criterion would imply that simulations that can estimate reliably the amount of dark mass 
within a couple of exponential scalelengths may be safely used for comparison with observations. 
Por bright spirals like the Milky Way this corresponds to radii of about 5-10 kpc, well outside 
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the region that may be compromised by numerical artifacts in the current generation of N-body 
experiments. 

In this paper we follow this proposal and compare the results of recent very high-resolution 
simulations of the formation of dark halos in the sACDM model with observational constraints 
on the total dark mass within spiral disks derived from the Tully-Fisher relation and from 
observations of Galactic dynamics. The numerical setup of the simulations is identical to that 
described by Navarro, Prenk & White (1997, hereafter NFW97) but the number of particles has 
been increased more than tenfold. As a result, each simulated sACDM halo has of order 250, 000 
dark matter particles within the virial radius and several thousands within radii comparable to 
the Sun's distance from the Galactic center (the "solar circle" Rg) so our numerical uncertainties 
are for all practical purposes negligible. 

We show below that comparison between these simulations and observational constraints 
reveals a severe inconsistency: sACDM halos have substantially more mass near the center than 
the maximum inferred from observations. Wc argue that this presents a serious challenge to 
the sACDM cosmogony and to many of its likely variants and that possible solutions involve 
fundamental revision of the basic premises or parameters of the model. 

2. Observational Constraints on Dark Mass in Individual Galaxies 

2.1. The Milky Way 

Kinematic observations of stars and gas in the Galaxy provide tight constraints on the total 
amount of dark matter within the solar circle, Ro- A direct estimate can be made by assuming 
that the halo is spherically symmetric, Mdark(?' < Ro) = ^davkiRo) Ro/G, where Vdark(-Ro) is the 
contribution of the dark halo to the circular velocity at Ro- This may be obtained by subtracting 
the disk contribution from the total circular velocity, V^^ark(^o) = ^c{Ro) — ^disk(-^o)- For the 
lAU-sanctioned values of Rq = 8.5 kpc and Vc{Ro) = 220 km s^^, and assuming that the disk 
potential is well approximated by an exponential disk with scalelength rdisk = 3.5 kpc and total 
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mass Mfjisk = 6 x (Binney & Tremaine 1987), we find 

Mdark(r- < Ro) = 5.2 X IQIOMq. (1) 

The uncertainty in this determination is hard to assess, although the evidence suggests that the 
mass in eq. 1 is actually an upper limit to the dark mass inside the solar circle. This is in good 
agreement with the recent Milky Way mass models of Dehnen & Binney (1998), who find that, 
within the slightly larger radius of 10 kpc the dark halo accounts for less than about 5-6 x W^^Mq, 
and perhaps as little as 2.5 x W^^Mq. 

The disk contribution to the circular velocity increases with our estimated value of -Mdisk, which 
in turn depends on: (i) the local density of the disk derived from the vertical kinematics of stars in 
the solar neighborhood (i.e., from "Oort limit" analysis, ~ 70Mqpc~^), (ii) on the exponential 
scalelength of the disk, and (iii) on the solar circle itself, through Mdisk oc r^j^j^ 6^°^'''*'=''. Recent 
reviews (Reid 1993, Sackett 1997) of available data suggest that the exponential scalelength 
assumed above may need to be revised downwards by up to 20%, but otherwise leave Sq and Rq 
largely unchanged]^ from the values assumed above. Such revision would increase the disk mass, 
leading to values of Mdark('^ < Ro) lower than quoted in eq. 1. Indeed, Sackett (1997) concludes 
that the Milky Way disk may very well be "maximal" once this revision is taken into account. 
Finally, we note that our procedure neglects the contribution of the Galaxy's bulge, lending further 
support to our interpretation of eq. 1 as an upper limit to the dark mass inside Ro- 

Figure 1 compares, as a function of halo mass, the dark mass estimate in eq. 1 with the 
results of simulations of several sACDM halos. Halo masses (M200) are measured inside the radius, 
^200) of a sphere of mean density 200 times the critical density for closure, and are typically 
characterized by the circular velocity at that radius, V200 = {GM2oo/r2oo)^^'^ = (W G Hq M2oo)^^^ ■ 
The reason for this choice is that at r2oo the circular orbit timescale is approximately equal to the 

^Note, however, that Oiling & Merrifield (1998), however, argue for Ro ^ 7.1 kpc upon analysis 
of the role of interstellar gas on the local values of Oort's constants. Since the circular velocity at 
Ro decreases correspondingly, this modification has little effect on our conclusions. 
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age of the universe; r2oo thus separates the "viriahzed" region of the halo from the region where 
mass shells are infalling into the system for the first time. The relevance of this definition stems 
from the latter property: only baryons inside r2oo may contribute to the baryonic mass of the 
galaxy, since those beyond this radius have yet to reach the center of the halo (White et al 1993). 

This property may be used to derive a firm lower limit to the mass of the halo that 
surrounds the Milky Way, corresponding to the case where the baryonic mass of the disk 
equals the total baryonic mass inside r2oo: i.e., Mdisk < -^dLIf ~ (^6/^o)-^200) implying 
that ^200 ^ [lOGHo (J^oM) Mdisk]^/^- For the disk mass adopted above and the cosmological 
parameters appropriate to sACDM discussed in §1, we find that 

F200 >130kms-^ (2) 

We emphasize that this is a strict lower limit to the mass of the halo that surrounds the Milky 
Way, since experiments show that typically not more than 80% of the baryons within r2oo are 
actually accreted into the central disk (Navarro &; White 1994, Navarro &; Steinmetz 1997). A 
similar, albeit more stringent, constraint may be derived by comparing the angular momentum of 
the Milky Way with those of sACDM halos. Halos with V200 < 150 km s~^ typically have specific 
angular momenta lower than that of the Milky Way disk (Mo, Mao & White 1998, Syer et al 
1999) and are therefore unlikely hosts of the Galaxy. Indeed, baryons typically lose a significant 
fraction of their angular momentum as they collapse to the disk (Navarro & White 1994, Navarro 
&; Steinmetz 1997) so we conclude that almost certainly the circular velocity of the Milky Way 
halo must exceed 130-150 km s~^. 

Figure 1 compares the constraints from eqs. 1 and 2 with the results of our sACDM numerical 
simulations. The comparison shows clearly a major discrepancy between the maximum dark matter 
inside Rg allowed by observations and the results of the numerical experiments. For example, 
sACDM halos with circular velocities similar to that of the Milky Way disk (V200 ~ VciRo) = 220 
km s^^) have about three times more dark mass inside the solar circle than inferred from 
observations. Even for the extreme case where the halo has the strict minimum circular velocity 
allowed by eq. 2, the simulations indicate an excess of more than 50% in the dark mass within Rg. 
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This serious discrepancy only worsens if we take into account that some extra dark material 
may be drawn inside Ro by the formation of the disk. A rough estimate of the magnitude of this 
correction can be made by assuming that the halo responds adiabatically to the assembly of the 
disk; the discrepancy then increases from 50% to almost 80% for the least massive halo allowed by 
eq. 2. We conclude that halos formed in the sACDM scenario are too centrally concentrated to be 
consistent with observations of the dynamics of the Galaxy. 

2.2. The TuUy-Fisher relation 

It is possible to extend the analysis of the previous subsection to a large fraction of disk 
galaxies by examining the tight correlation between the total luminosity of galaxy disks and the 
rotation speed of their gas and stars (the Tully- Fisher relation, Tully Sz Fisher 1977). Provided 
that disk mass-to-light ratios and exponential scalelengths can be estimated reliably, it is possible 
to evaluate the disk contribution to the circular velocity and to apply the same analysis of the 
previous subsection to derive constraints on the total dark mass contained within the optical 
radius of the galaxy. 

We choose to carry out the analysis at 2.2 exponential scalelengths (2.2rdisk) from the 
center, since the contribution of exponential disks to the circular velocity peaks there and it is at 
that radius that TuUy-Fisher velocities are typically measured (Courteau 1997). (We note that 
Ro ~ 2.4rdisk in the case of the Milky Way, so this choice of radius is similar to that adopted in 
§2.1.) Constraints on dark masses inside 2.2r(jisk depend sensitively on estimates of the total mass 
associated with galaxy disks. Because of the lack of an Oort-limit analog in external galaxies, we 
resort to broad-band colors as estimators of the mass-to-light ratio of the stellar disk. Late type 
spirals such as those that make up the majority of galaxies in Tully-Fisher samples have B — R 
colors in the range (0.2, 1.0) (Courteau 1999) which, for a galaxy that has been steadily forming 
stars for ~ 13 Gyrs, imply /-band mass-to-light ratios of order (Mjisk/L/) ~ 2 it 1(M/L/)q. This 
estimate assumes a Salpeter initial mass function and exponential star formation histories with 
timescales which vary from tsf ~ 1 Gyr to star formation rates that are constant over the age of 
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the universe. Similar values are obtained using the GISSEL96 models of Bruzual & Chariot (1996) 
and the PEGASE models of Fioc & Rocca-Volmerange (1997). Uncertainties in this estimate of 
{Mdisk/ Li) are hard to assess, but are unlikely to be larger than about a factor of two, which is 
the value of the "error" we assume here. ^ 

With this caveat, we compute the rotation speeds (14ot = ^4(2.2 rrfj^fc)) of hypothetical 
exponential disk galaxies with Mdisk/Lj = 2(M/Lj)q assembled at the center of the simulated 
dark halos shown in Figure 1. Disk radii are chosen so as to satisfy the empirical relation, 
^-disk ~ 3(I4ot/200 km s~i)/i-i kpc (Courteau 1999, Mo, Mao & White, 1998, Navarro 1999). 
Since disk radii depend on Vj-ot, an iterative procedure is needed, which we implement as follows. 
Given a halo of circular velocity V200 we assign to it a disk with exponential scalelength derived 
assuming that Vrot = ^00 • The simulated halo structure and the disk potential are then used to 
compute a new Vj^ot estimate, taking into account the "adiabatic" response of the halo to the disk 
potential (see Navarro, Frenk & White 1996, hereafter NFW96, for details). This new velocity 
estimate is used to recompute rjisk and the procedure is then iterated until convergence. 

The solid-line curves in Figure 2 illustrate the result of applying this procedure to three 
representative dark halos of different mass (or luminosity, since we assume a constant mass-to-light 
ratio) , as a function of the total mass adopted for the disk of the galaxy. In each case we vary the 
total disk mass from zero to M^^^, the maximum value compatible with the baryonic content of 
the halo. As the disk mass increases, each hypothetical galaxy moves from left to right across the 
plot. When the disk mass becomes comparable to the dark mass inside 2.2r(iisk the curve inches 
upwards and becomes essentially parallel to the observed Tully-Fisher relation. The increased 

^We note that the same procedure gives F-band mass-to-light ratios of about 2{M/ Lv)q for 
galaxies in the same color range, more than a factor of 2 lower than that derived for the Milky Way 
disk from the Oort limit (Binney &: Tremaine 1987). This reflects the well known Galactic "disk 
dark matter" content. Assuming that Tully-Fisher galaxies have similar amounts of (presumably 
baryonic) "disk dark matter" would make all the arguments we develop here even stronger. 
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potential due to the disk mass and the extra dark mass drawn inside 2.2r(jisk under the assumption 
that the halo responds adiabatically to the assembly of the disk contribute similarly to the overall 
increase in rotation speed over and above that of the original dark halo. The rightmost point of 
each curve is reached when the mass of the disk equals M^^^. 

It is clear from Figure 2 that, even under the extreme assumption that galaxies contain all 
available baryons in each halo, simulated disks are almost two magnitudes fainter than observed. 
Increasing the baryonic mass of a halo has virtually no effect on this conclusion, since in this case 
model galaxies would just move further along paths approximately parallel to the Tully-Fisher 
relation, as shown in Figure 2. Disk galaxies assembled inside sACDM halos therefore cannot 
match the observed Tully-Fisher relation, unless one or more of the assumptions in our procedure 
are grossly in error. 

Perhaps the most uncertain step in our argument is the stellar mass-to-light ratio adopted 
for the analysis. The horizontal "error bar" shown on the starred symbols in Figure 2 indicates 
the effect on our results of varying the /-band mass-to-light ratio by a factor of two from our 
fiducial value of 2 in solar units. This is not enough to restore agreement with observations, which 
would require (Mdisk/-^/) ~ 0.4, a value much too low to be consistent with standard population 
synthesis models. The vertical "error bars" illustrate the effect of varying the "concentration" of 
each halo by a factor of two.^ Even with this large variation in halo structure, our hypothetical 
disks fail to reproduce the observations. 

A second uncertainty comes from the "adiabatic contraction" correction applied to dark 
halos in order to mimic the halo response to the disk assembly, and one may wonder whether our 
"adiabatic" contraction assumption is at all correct. This possibility may be checked, however, 
through direct numerical simulation. We have therefore included gas in the simulations and have 

^The NFW "concentration parameter" is defined as c = r2oo/rs, where the scale radius, 
r<j is one of the parameters of the density profile model proposed by NFW96 and NFW97, 
p{r) oc (r/rs)^^(l -|- r/r^)"^, see those papers for details. 
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evolved them again using with the simpHfied treatment of radiative coohng and star formation 
described in detail in Steinmetz Sz Navarro (1999). The results are shown as open squares in 
Figure 2, and are in good agreement with the results of the simple modeling proposed above. 

As in the case of the Milky Way discussed in §2.1, the problem can be traced to the large 
central concentration of sACDM halos. Indeed, one way to solve the problem would be to 
reduce substantially (by a factor of two or three) the dark mass in the innermost few kpc of 
galaxy-sized dark halos. This would significantly reduce T4ot) bridging the gap between model 
and observations. In terms of the halo density profile model proposed by NFW96 and NFW97, 
this would be equivalent to reducing the "concentration parameter" , c, by a factor of about five. 
This is the same conclusion reached by Navarro (1999), who advocated that halos with very low 
"concentration parameters", c<3, and stellar mass-to-light ratios as low as 0.5 /i(M/L/)q in 
low surface brightness galaxies were required in order to match the shapes of disk galaxy rotation 
curves. We emphasize, however, that the present conclusion is independent of assumptions about 
the detailed shape of the dark matter density profile, and depends largely on the total dark mass 
on scales of individual galaxies. 



3. Discussion and Conclusions 

The analysis of the previous subsection demonstrates conclusively that the central mass 
concentration of sACDM halos is inconsistent with observations of the dynamics of spiral galaxies. 
Does this argument rule out all models based on the CDM-|-inflation paradigm or are there 
plausible modifications to the sACDM parameters that may bring the model into agreement with 
observations? 

In principle, there are parameter choices that may help make ACDM models consistent with 
observations of the internal dynamics of galaxies, but the magnitude of the modifications required 
are, however, uncomfortably large, and come at the expense of other major successes of the model. 
As discussed by NFW97, the central concentration of dark halos is directly proportional to the 
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mean matter density of the universe at a suitably defined collapse time. For halos of fixed mass, 
these collapse times may depend sensitively on the values adopted for the cosmological parameters. 
However, as discussed by Navarro (1999, see his Figure 8), the combination of parameters needed 
to reproduce the present-day abundance of galaxy clusters is such that the characteristic densities 
of galaxy-sized dark halos is approximately independent of JIq and of the value of Hq. Indeed, 
halos in the former "standard" CDM model (JIq = 1, h = 0.5, as ~ 0.6) have very similar 
concentrations as the sACDM halos we discuss here since both models are normalized to match 
the abundance of massive clusters at z = 0. This implies that all CDM models that match the 
abundance of clusters are likely to have difficulty reproducing the dynamics of spiral galaxies. 

Even if cluster normalization is dropped from the list of relevant constraints and only the 
COBE measurements are used to normalize the power spectrum of initial density fluctuations, 
reducing substantially the central concentrations of dark halos still implies uncomfortable 
parameter choices. For example, if JIq is taken as a free parameter, reducing the dark mass of 
a V200 ^ 200 km s"^ halo inside Rg ~ 8.5 kpc by a factor of three relative to sACDM requires 
ft() ~ 0.05, a very low number that would call into question the need for a dominant non-baryonic 
dark matter component in the Universe. Similarly, if rather than Qq, Hq is allowed to vary, values 
as low as 30 km s~^ Mpc~^ are required to obtain the desired effect, in gross disagreement with 
current observational estimates. 

If wc arc to preserve the successes of the sACDM model, what seems to be required is 
a substantial change in the shape of the power spectrum relative to that predicted by the 
CDM+inflation paradigm. A "designer" power spectrum that would go some way towards 
reconciling the model with observations would suppress power on galactic and subgalactic scales 
while keeping the large scale properties of the model virtually unchanged, as envisioned, for 
example, in models that introduce a "tilt" in the primordial power spectrum relative to the 
standard Harrison-Zel'dovich value. This would in principle allow galaxy-sized dark halos to 
collapse later and thus become less centrally concentrated, although the magnitude of the tilt 
required is still unclear. Other alternative modifications that may in principle reproduce the 
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desired trend would involve the existence of a sizeable "hot dark matter" component (although 
much higher than derived from current estimates of the neutrino mass, see, e.g., Primack & Gross 
1998), or possibly even dark matter candidates that may annihilate without trace in dense regions 
such as the centers of galaxies. One problem that afflicts all these proposed modifications is that 
they may hinder the formation of massive galaxies at high redshift, at odds with the mounting 
evidence that such galaxies are fairly common at z>3 (see, e.g., Steidel et al 1998). A thorough 
investigation of these possibilities will be needed in order to assess just how radically the new 
"standard" model of structure formation must be revised in order to bring spiral galaxies into the 
realm of observations that are consistent with the current paradigm of structure formation. 

This work has been supported by the National Aeronautics and Space Administration under 
NASA grant NAG 5-7151 and by NSERC Research Grant 203263-98. MS acknowledges useful 
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Simon White, Carlos Frenk, Mike Balogh, and Mike Hudson, as well as the hospitality of the 
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Fig. 1. — Dark mass enclosed within a radius Ro = 8.5 kpc, the Sun's distance from the center of the 
Milky Way, versus the circular velocities of sACDM halos. The shaded region highlights the allowed 
parameters of the dark halo surrounding the Milky Way, as derived from observations of Galactic 
dynamics and by assuming that the disk mass cannot exceed the total baryonic content of the 
halo. The filled circles show the loci of sACDM halos as determined from high-resolution N-body 
simulations. The solid line is the circular velocity dependence of the dark mass expected inside Rq 
for halos that follow the density profile proposed by NFW96 and NFW97. The circular velocity 
dependence of the NFW "concentration" parameter of the simulated halos is well approximated on 
these scales hy c ^ 20 (V20o/100kms^^)^"'^/^ (dotted line). This is slightly higher than predicted 
by the approximate formula proposed by NFW97 but consistent with their published results. 
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Fig. 2. — The /-band TuUy Fisher relation compared with the loci of hypothetical exponential 
disk galaxies assumed to assemble at the center of three representative sACDM halos. Dots are 
a compilation of the data by Giovanelli et al. (1997), Mathewson, Ford Sz Buchhorn (1992) and 
Han & Mould (1992). The solid line is the best fit to the data advocated by Giovanelli et al. The 
hypothetical galaxies have radii consistent with observations and move from left to right along each 
curve (labeled by the circular velocity of the halo at the virial radius) as the disk mass increases, 
under the assumption of a constant stellar mass-to- light ratio, M/Lj = 2 in solar units. The starred 
symbols correspond to the maximum disk mass, M^^^, allowed by the universal baryon fraction 
of the sACDM model. Open squares are N-body gasdynamical simulations of the formation of 
galaxies within these halos. Error bars correspond to two different choices of IMF, as discussed in 



